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Abstract. We report high angular resolution (HPBW~ 0.6" x 0.5" at 1.3mm) observations of the Class 
intermediate-mass (IM) protostar NGC 7129-FIRS 2 using the Plateau de Bure Interferometer. Our observa- 
tions show the existence of an intense unresolved source in the continuum at 1.3mm and 3mm at the position 
of the Class object. In addition, compact CH3CN emission is detected at this position. Fhe high rotational 
temperature derived from the CH3CN lines (T rot «50 K) as well as the enhanced CH3CN fractional abundance 
' (X(CH3CN)~7.0 10~ 9 ) show the existence of a hot core in this IM young stellar object. Fhis is, up to our knowl- 

■^J- , edge, the first IM hot core detected so far. Interferometric maps of the region in the CH3OH 5kk/^4 kk i , and D2CO 

' 4o4— >3o3 lines are also presented in this paper. Fhe methanol emission presents two condensations, one associated 
with the hot core which is very intense in the high upper state energy lines (E„>100 K) and other associated 
with the bipolar outflow which dominates the emission in the low excitation lines. Enhanced CH3OH abundances 
(XfCIFjOH) ~ 3 10~ 8 - a few 10 -7 ) are measured in both components. While intense D2CO 4o4^3o3 emission 
is detected towards the hot core, the N2D + 3— >2 line has not been detected in our interferometric observations. 
Fhe different behaviors of D2CO and N2D + emissions suggest different formation mechanisms for the two species 
|' and different deuteration processes for H2CO and N2H + (surface and gas-phase chemistry, respectively). Finally, 

1~ \ the spectrum of the large bandwidth correlator show a forest of lines at the hot core position reavealing that 

this object is extraordinarily rich in complex molecules. To have a deeper insight into the chemistry of complex 
molecules, we have compared the fractional abundances of the complex O- and N- bearing species in FIRS 2 with 
those in hot corinos and massive hot cores. Within the large uncertainty involved in fractional abundance esti- 
mates towards hot cores, we do not detect any variation of the relative abundances of O- and N-bearing molecules 
([CH3CN]/[CH3 0H]) with the hot core luminosity. However, the O-bearing species H2CO and HCOOH seem to 
be more abundant in low and intermediate mass stars than in massive star forming regions. We propose that this 
could be the consequence of a different grain mantle composition in low and massive star forming regions. 

Key words. Stars: formation - individual: NGC7129 - FIRS 2-ISM: abundances - clouds - individual: NGC 7129 



1. Introduction 

Hot cores are compact objects near or around protostars 
characterized by warm temperatures (T\.>100 K) and 
high densities (n>10 6 cm -3 ). These regions are also char- 
acterized by a very rich chemistry in complex molecules 
(CH3OH, CH 3 CN, CH3OCHO, CH 3 OCH3,C 2 H 5 CN ...). 
Hot cores are thought to be associated with high- 
mass protostars (M>8 M©) and to represent an im- 
portant phase in their evolution toward ultracompact 
and compact HII regions. Recently, regions character- 
ized by warm temperatures and high densities have also 
been detected in two low-mass protostars IRAS 16293- 
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2422 <|Ceccarelli et al. 2000t ICazaux et al. 2003h and 
NGC 1333 IRAS 4A l|Bottinelli et al. 2004fh Complex 
molecules typical of hot cores (e.g. HCOOH, CH3OCHO, 
CH 3 CN, C2H5CN) have also been detected in these ob- 
jects. However, the amount of warm material involved, as 
well as the chemistry, are different in the two classes of 
objects. For this reason, the warm regions in the inner en- 
velope of low mass protostars are usually referred to as 
"hot corinos". 

The formation of complex molecules in hot cores and 
corinos is poorly understood. In the standard scheme, 
neutral molecules (CO, CS..) are frozen onto dust grains 
during the cold pre-stellar phase. If the dust tempera- 
ture is sufficiently low during this phase, surface hydro- 
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genation of CO leads to the formation of solid H 2 CO 
and CH 3 OH (e.g. Brown ct al. 1988; Charnley et al. 
1992; Caselli et al. 1993). Once the star starts heat- 
ing the grain surfaces, these molecules (called "parent" 
species) evaporate, enlarging their abundances in the gas 
phase. Because of the high temperature and density of hot 
cores and corinos, these molecules undergo fast neutral- 
neutral and ion-neutral reactions producing at early stages 
(~ 10 4 yr) a second generation of complex O-bearing 
species called "daughter" molecules, such as methyl for- 
mate, CH3OCHO. However, Horn et al. (2004) recently 
found that the gas-phase reaction sequence between pro- 
tonated methanol and formaldehyde, crucial for the gas 
phase formation of methyl formate, does not proceed in 
their laboratory experiments. This suggests that surface 
chemistry is probably playing a key role also in the for- 
mation of this species. 

Regarding complex N-bearing species (e.g. CH3CN, 
CH3CH2CN), they are observed in hot cores and hot cori- 
nos. In the chemical scheme of Rodgers & Charnley (2001), 
complex N-bearing species are thought to be formed in the 
gas phase about 10 5 yr after the grain mantle evaporation, 
only if the gas temperature is sufficiently large (> 300 K). 
On the other hand, Caselli et al. (1993) found that large 
abundances of ethyl cyanide (X(CH 3 CH 2 CN) - 10~ 7 , 
w.r.t. H 2 ) can be obtained on grain surfaces, if the dust 
temperature during the protostellar accreting is around 
40 K. In the case of lower temperatures (10-20 K), the 
CH3CH2CN can still be formed on the surface, but with 
reduced abundances (~ 10~ 9 ). Thus, Caselli et al. (1993) 
suggests that saturated (or H-rich) complex N-bearing 
species (such as ethyl cyanide) are "parent" species and 
that they should be observable as soon as the grain man- 
tles are released in the gas phase upon stellar heating of 
the dust. Later on, "H-poor" N-bearing molecules (e.g. 
vinyl cyanide, CH2CHCN) can be formed in the gas phase 
from the destruction of CH3CH2CN. 

Indeed, if complex species in general are not formed 
on grain surfaces, there are important problems for un- 
derstanding the chemical composition of hot corinos: the 
timcscale necessary to convert "parent" molecules into 
complex "daughter" molecules is much longer than the 
transit time of the gas in hot corinos (a few hundred 
years). Therefore, it seems that the chemistry of complex 
molecules must begin on grain surfaces (at least in hot 
corinos). However, chemists are still far from a unique in- 
terpretation of the hot core (and corino) chemistry. More 
observational work is needed. In particular, given the large 
sensitivity of surface and gas-phase chemistry on dust 
temperature and gas density, it will be extremely impor- 
tant to measure possible variations of complex molecule 
abundances with the physical characteristics of hot cores, 
which of course depend on the mass and luminosity of the 
associated protostar. 

In this paper, we present an interferometric study of 
the hot core associated with the Class IM protostar 
NGC 7129-FIRS 2. With a luminosity - 500 L Q and a 
mass ^5 M , FIRS 2 is very likely the youngest IM ob- 



ject known at present (Fuente et al. 1998,2002). An ener- 
getic bipolar molecular outflow with a quadrupolar mor- 
phology is associated with it (Fuente et al. 2001, hereafter 
Paper I). Interferometric observations in the continuum at 
1mm and the 12 CO 2— >1 line shows that the quadrupolar 
morphology of the outflow is due to the superposition of 
two bipolar molecular outflows FIRS 2-out 1 and FIRS 2- 
out 2. FIRS 2-out 1 is associated with the Class proto- 
star detected as an intense millimeter clump in the contin- 
uum image (MM1 in the nomenclature of Paper I) while 
FIRS 2-out 2 is associated with a more evolved infrared 
star (FIRS 2 - IR) undetected at millimeter wavelengths. 

A quite complete chemical study of this Class IM 
source has been carried out by Fuente et al. (2005) (here- 
after Paper II) using the 30m IRAM telescope. They de- 
tected warm CH 3 CN ( Tfe>63 K) towards the central po- 
sition which constitutes a strong evidence for the existence 
of a hot core in this Class IM object. However, the lim- 
ited angular resolution of these observations make it very 
difficult to discern between the hot core emission and those 
of the outflow and/or the warm envelope. The new results 
presented in this paper confirm the existence of the hot 
core associated with FIRS 2 and give a first glance of the 
chemistry of this object. Up to our knowledge, FIRS 2 is 
the first IM hot core detected thus far. The intermediate 
kinetic temperature and mass of the FIRS 2 hot core are 
expected to produce a differentiated chemistry and to fur- 
nish a link between the low mass and high mass regimes. 



2. Observations 

The observations presented in this paper correspond to 
two different observational projects carried out with the 
Plateau dc Burc interferometer (PdBI). The main set of 
data was observed on March 08, 2003. The observations 
were carried out with excellent weather conditions and 
counts 5 hr of integration time (on-source) with 6-elemcnt 
array in the A configuration which provides the highest 
angular resolution. The spectral correlator setup was ad- 
justed to observe the CH3CN 5^4 line with two contigu- 
ous 20 MHz units, N 2 D+ 3^2 and D 2 CO 4 04 ^3 3 with 
two separated 20 MHz units and the remaining units were 
configured for maximum continuum sensitivity (see Table 
1). These observations show on average system tempera- 
tures of 140 K (3mm) and 250 K (1mm), a mean atmo- 
spheric water vapour content of l-2mm, and were made 
under seeing conditions of 0.3 arcsec. The RF calibrator 
was 0420-0414, the phase calibrators ordered by right as- 
cension were 1928+738, 2146+608 and 2309+454, and the 
flux calibrator was MWC 349. The precision in the flux 
densities is better than 20% at 1mm and better than 10% 
at 3mm. Images have been created in natural weighting. 

In addition, we present the CH 3 OH J=5fc— >4fc im- 
age observed within a previous project in November and 
December, 1998. The data corresponding to this project 
were partially published in Paper I. Details about the ob- 
servations are given in Table 1 and Paper I. 
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Fig. 1. Intcrfcromctric maps of the Class protostar 
NGC 7129-FIRS 2 in the continuum at 1mm and 3mm, 
and the following molecular lines: CH3CN 5^4, D2CO 
4o4— >3o3, CH3OH 51.4^4^4. The panel marked with 
"CH3OH (high)" show the integrated intensity map of 
the lines 5 3) i— >4 3) i, 5 3 . 2 ^ 3 . 2 ,5 2 ,2^4:2,2 and 5 3 ,3->4 3i 3. 
Contour levels are: a) 10 to 60 by 5 mJy/beam; b) 20 to 
160 by 20 mJy/beam ;c) 0.6 to 2.7 by 0.3 Jy/beam km s _1 ; 
d) 0.3 to 0.6 by 0.1 Jy/beam km s" 1 ; e) 1 to 4 by 0.5 
Jy/beam km s _1 ; f) 1 to 7.5 by 0.5 Jy/beam km s -1 . 
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Fig. 2. Continuum visibilities at 92 GHz and 230 GHz vs 
projected baselines in units of wavelength. In the bottom 
panel we show together the 92 GHz and the 230 GHz visi- 
bilities. The 92 GHz visibilities have been scaled by a fac- 
tor of 10.5 to match the 230 GHz ones. The perfect match 
between the 92 GHz and 230 GHz visibilities suggests that 
both emissions arise in the same region as expected in the 
case of dust thermal emission. 



3. Continuum data 

The continuum images at 3mm and 1.3mm are shown 
in Fig. 1. The angular resolution of these images is 
3 times better than that of the images previously 
published in Paper I. Both images show an intense 
and compact source located at RA(2000)=21:43:01.684, 
Dec(2000)=66:03:23.619, the position named MM1 in 
Paper I. There is no indication of the presence of the sec- 
ond and more extended source MM2 which is very likely 
resolved out in the high angular resolution image. 

We have modeled the visibilities to have a deeper in- 
sight into the source structure. In Fig. 2 we show the 
92 GHz and 230 GHz visibilities vs projected baselines 
in units of wavelength. The 92 GHz visibilities scaled by 
a factor of 10.5 match perfectly the 230 GHz visibilities. 
This suggests that both emissions are arising exactly in 
the same region. The spectral index is a=2.56. The fact 



that the 3mm and 1.3mm emission arises in the same re- 
gion and the spectral index larger than 2 is consistent with 
dust thermal emission. We have used the high-angular 
resolution 1.3mm image to model the continuum emis- 
sion. Several models (elliptical Gaussian, disk, F„=^~ 2 , 
Fv=v~ 3 , and an elliptical Gaussian+point source) have 
been used to fit the visibilities. The best fit has been ob- 
tained with the elliptical Gaussian+point source model. 
The HPFW of the elliptical Gaussian is -650x900 AU 
and the flux is 0.43 Jy. The point source has a flux of 
0.13 Jy. The total flux of the compact 1.3 mm compo- 
nent is ~0.56 Jy. This implies a total (dust +gas) mass 
of 2 M Q assuming a standard dust temperature of 100 K 
and a dust emissivity fc„=0.015(1300/A (^m)) cm 2 g _1 . 
Taking into account the uncertainty in the value of the 
dust temperature which could range between 50 K and 
300 K this mass is accurate within a factor of 4. 
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Table 1. Observations with the Plateau de Burc 
Interferometer (PdBI) 





Conf. 


Beam 


date 


Continuum 3mm 


A 


1.56" x 1.20" 


March 03 


Continuum 1mm 


A 


0.63" x 0.46" 


March 03 


CH 3 CN 5 fc ^4 fc 


A 


1.56" x 1.20" 


March 03 


N 2 D+ 3^2 


A 


0.63" x 0.46" 


March 03 


D 2 CO 4 4^3o3 


A 


0.63" x 0.46" 


March 03 


CH3OH 5 M ,->4 fc ,v 


CD 


1.51" x 1.42" 


Nov, Dec 98 



We speculate about the possibility that the point 
source is an accretion disk. In this case, the mass of the 
disk would be ^0.3-0.8 M Q . This value is 30 times larger 
than that found by Fuente et al. (2003) in the Herbig Be 
star R Mon, but it is similar to that found by Rodriguez et 
al. (2005) in the Class low-mass protostar IRAS 16293- 
2422B. Thus, if our assumption is confirmed, the differ- 
ence in the disk occurrence and masses between IM and 
low mass stars would be related to a short timescale for 
the disk dissipation instead of to differences in the first 
stages of the star formation. This interpretation is in line 
with Fuente et al. (2003). 

4. Molecular line images 

4.1. CH 3 CN 

The interferometric image of the Class protostar 
NGC 7129-FIRS 2 in the CH 3 CN 5^4 line shows a com- 
pact source at the same position and with the same size 
as the 1mm continuum source (Fig. 1). The total agree- 
ment in position and size between the continuum and the 
CH3CN emissions points out to the existence of a well- 
differentiated compact component in the protostellar en- 
velope. 

Because of the rotational structure of CH3CN, one can 
observe several lines at different energies very close in fre- 
quency. In Fig. 3 we show the interferometric spectrum of 
the CH3CN 5^4 line compared to the single-dish spec- 
trum reported in Paper II. The interferometer has recov- 
ered 100% of the flux of the highest upper state energy 
component (E u =124 K) (see Fig. 3). However, only 37% 
of the flux of the lowest energy component (E u =13 K) has 
been recovered (see Fig. 3). This suggests that while the 
emission of the low energy CH3CN lines mainly arises in 
the cooler and more extended envelope, the emission of 
the high energy lines entirely arises in the hot core. In ad- 
dition to the CH3CN K-ladder we have also detected the 
CH3 13 CN 5o~^4o line which allows us to have an estimate 
of the opacity of the CH3CN line. 

The CH3CN column density has been estimated using 
the rotational diagram technique. The low energy lines are 
expected to be optically thick in the hot core. This is con- 
firmed by our interferometric observations. In particular, 
we measure I(CH 3 CN 5 ^4 )/I(CH 3 13 CN 5 ^4 )~2.4 
which implies an opacity ~25 for the main isotope line 



assuming 12 C/ 13 C=70. We have corrected by the opacity 
effect and derived T rot '^54 K and a beam averaged col- 
umn density N(CH 3 CN)~1.1 10 16 cm" 2 for the hot core 
component using the 5q-^4o and 5s^4s lines. The de- 
rived rotation temperature is in agreement, within the 
uncertainties, with that derived from single dish data 
l|Fuente et a l. 2005). Since the low energy lines are opti- 
cally thick in the hot core, we think that this agreement 
is fortuitous and due to the addition of two uncertain- 
ties which shift the estimated column density in opposite 
directions, the contribution of extended emission in the 
single-dish beam and the opacity effect. From our inter- 
ferometric observations, we obtain X(CH3CN)^7.0 10 -9 
in the 2 M Q hot core in FIRS 2 assuming a source size of 
0.72" xO. 52" (see Table 5). This fractional abundance is 
similar to that measured in massive hot cores and slightly 
larger than those derived in hot corinos. 

For comparison, we have estimated the CH3CN abun- 
dance in the cooler and extended envelope. For this aim, 
we have subtracted the emission of the compact compo- 
nent from the single-dish spectrum and analyze the re- 
sult using a rotational diagram. The opacity effect are 
expected to be less important in the low density envelope 
than in the hot core. We derive T ro t=12 K and a beam 
averaged column density of 5.8 10 11 cm" 2 (HPBW=25"). 
Using the H2 column density derived from single-dish con- 
tinuum observations by Fuente et al. (2001), we estimate a 
CH3CN fractional abundance of ~3 10~ n in the extended 
component, i.e. about 3 orders of magnitude lower than 
in the compact component. 

Summarizing, our CH3CN interferometric data show 
the existence of a compact (900x650 AU) source charac- 
terized by a high kinetic gas temperature (>50 K) and en- 
hanced CH 3 CN abundance (X(CH 3 CN)~7.0 10~ 9 ) in the 
Class IM protostar FIRS 2. This constitutes a definite 
proof of the existence of a hot core in this IM protostar 
which, up to our knowledge, it is the first IM hot core 
detected so far. 

4.2. CH 3 OH 

We have observed the CH3OH ik,k'=5k,k>^^k,k> lines 
towards FIRS 2 using the 2.5 MHz wide band corre- 
lator which allows us to observe most of the CH3OH 
Jfe,fe'=5fe.fc'^'4fe ) fe' components simultaneously. In Fig. 1 
we show the total integrated intensity map of the CH3OH 
5^4— >4i.4 line which is that observed with the lowest up- 
per stage energy (E„=35 K). The emission of this com- 
ponent presents two maxima, the first one is spatially co- 
incident with the hot core while the second is located at 
the position (—2", — 1.5"). This position lies on the axis of 
the outflow FIRS 2-out 1. In Fig 1 we also show the inte- 
grated intensity map of the 5 3i i^4 3! i, 5 3! 2— >4 3! 2,52,2-^4 2: 2 
and 53 t 3— >43_3 lines [panel indicated as CH3OH (high) in 
Fig. 1]. These lines are partially overlapped in our spec- 
trum and all of them have upper state energies >50 K. 
In this case, the emission is mainly arising in the hot core 
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5 



Position" 



Integrated intensity 6 

(Jy) 



Size c 
(arcseconds) 



Continuum 3mm 
Continuum 1mm 

CH 3 CN 5 k ^4 k 
D 2 CO 4o4^3o l3 



21:43:01.7 66:03:23.6 

21:43:01.7 66:03:23.7 

21:43:01.7 66:03:23.6 

21:43:01.7 66:03:23.6 

21:43:01.7 66:03:23.6 



0.0565(0.0005) 
0.4278(0.0047) 
0.1276(0.0046) 
0.0734(0.0060) 
0.76(0.6) 



0.66(0.02)x0.54(0.02) 
0.72(0.01) x0.52(0.01) 
Point source (< 0.3") 
0.63(0.22) x0.46(0.19) 
0.58(0.24) x0.40(0.24) 



a The absolute positional precision is < 0.3". 

6 Systematic calibration errors are not included in the error budget associated to the integrated intensity 
c Deconvolved by the synthesized interferometric beam. 




Fig. 3. Comparison between the single-dish (thin histograms or red) and interferometric (thick histograms or black) 
spectra towards the hot core position. The interferometric spectrum has been scaled for an easier comparison. Note 
that we have recovered all the flux for the CH3CN line with the highest upper state energy. 



with a very weak contribution of the second maximum. 
Thus, the morphology of the CH3OH emission shows the 
existence of two components, a compact one which is spa- 
tially coincident with the hot core (hereafter hot core com- 
ponent) and an extended one which is very likely associ- 
ated with the outflow FIRS 2-out 1 (hereafter outflow 
component) . 



In Fig. 3 we compare the interferometric spectra for 
the hot core and outflow components with the single- 
dish spectrum at the central position. The two compo- 
nents present differences in their kinematics and excita- 
tion conditions. The emission of the CH 3 OH lines towards 
the hot core is centered at a velocity — 10±1 km s _1 , 
while the emission of the outflow component is centered at 
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-6±1 km s" 1 . The 12 C0 2->l and SiO 2->l maps of the 
outflow FIRS 2- out 1 reported in Paper II show that this 
velocity, -6±1 km s , is characteristic of the high velocity 
bullet Rl, reinforcing our interpretation of this emission 
as associated with the bipolar outflow. 

Different excitation conditions characterize the hot 
core and outflow components. This is clearly seen when 
one compares the relative intensities of the high and low 
energy lines in the methanol spectra towards the two stud- 
ied positions. All the CH3OH lines have similar inten- 
sities towards the hot core. However, high energy lines 
are a factor of ~5 weaker than the low energy ones to- 
wards the outflow position (see Fig. 3). We have derived 
the rotation temperature and the beam averaged CH 3 OH 
column density for the hot core using the rotational dia- 
gram technique (see Fig. 4). We obtain T rot =809 K and 
N(CH 3 OH)=1.5 10 17 cm- 2 . This temperature is quite 
high compared with that obtained from the CH3CN lines. 
Similarly to the case of CH3CN, the low energy lines are 
very likely optically thick and the derived rotation tem- 
perature is an upper limit to the actual one. To have an 
estimate of the uncertainty due to the unknown rotation 
temperature in the CH 3 OH column density, we have cal- 
culated the CH3OH column density from the flux of the 
54,3^4,3 line (E u =104 K) assuming T rot =50 K. We ob- 
tain N(CH 3 OH)=4.9 10 16 cm" 2 . Thus, we consider that 
our beam averaged column density estimate is accurate 
within a factor of 3. Assuming T rot =50 K, the derived 
CH3OH fractional abundance is 3.0 1CP 8 in the hot core. 
This abundance is similar to that measured in the proto- 
typical massive hot core OMC1. 

For comparison, we have estimated the rotation tem- 
perature and methanol column density for the outflow 
condensation (see Fig. 4). We have obtained T rot =19 K 
and N(CH 3 OH)=2.9 10 15 cm" 2 . The low rotational tem- 
perature of this component supports our interpretation 
of being associated with the outflow (postshocked mate- 
rial) instead of with the hot core. We have not detected 
a millimeter continuum counterpart for this condensa- 
tion in our interferometric 1mm map. In principle, this 
could be a sensitivity problem. From our interferometric 
image and assuming Td„ st =100 K we derive a 3xcr up- 
per limit of ~0.01 M Q for the mass of the condensation 
associated with the methanol clump. This would imply 
X(CH 3 OH)>3.8 10~ 7 in this clump. We can alternatively 
think that this clump is not the result of a hydrogen den- 
sity enhancement at this position. For example, it could 
be produced by a local enhancement of the CH3OH abun- 
dance because of the shocks associated with the bipolar 
outflow. To estimate the methanol abundance in this case, 
we calculate the H 2 column density from the single-dish 
13 CO observations reported in Paper II. At the velocities 
of the bullet Rl, the 13 CO column density averaged in 
the 30m beam is ^3 10 15 cm~ 2 . If we assume a uniform 
distribution of the molecular gas in the beam, we obtain 
that the CH3OH abundance is similar to that of 13 CO for 
which we assume a standard value of ^10~ 6 . Obviously, 
this is an upper limit to the actual abundance since we 
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Fig. 4. Rotational diagram of CH3OH towards the hot 
core and the outflow condensation. Note that the rotation 
temperature towards the hot core is much larger than to- 
wards the outflow condensation. 

expect some dumpiness. Thus, the CH3OH abundance in 
the outflow condensation could be larger than in the hot 
core and reach values close to 10 ~ 6 . This large value of 
the methanol abundance is consistent with those found 
in the molecular bullets associated with low-mass bipolar 
outflows (see e.g. Bachiller & Perez-Gutierrez 1997). 

In addition to the CH3OH lines, the spectrum towards 
the hot core show several lines that are weaker towards 
the outflow condensation. We have identified these lines 
as belonging to the complex molecules CH3OCH3 and 
C2H5CN. We have also a tentative detection of the sul- 
phuretted species 34 S02- The detection of complex satu- 
rated molecules corroborates the existence of a hot core 
in this source and reveals the rich and complex chemistry 
associated with it. 



4.3. D 2 CO 

The interferometric map of D2CO shows a compact source 
located at the position of the hot core (see Fig. 1). The size 
of the D2CO emitting region is similar to that of CH3CN. 
However, the linewidth of the D 2 CO line is ~11 km s _1 , 
larger than those of the CH3CN and CH3OH lines which 
are typically ~6-8 km s . Assuming the standard inside- 
outside collapse model, the largest linewidth of the D 2 CO 
line is consistent with the emission arising in an inner 
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region of the protostcllar envelope. In a simple calculation, 
an accretion velocity of ~5 km s" 1 corresponds to a radius 
of ^130 AU assuming a dynamical age of ^10 5 yr and 
dM/dt= 5 10~ 5 M Q yr -1 . A similar radius was found 
by Maret et al. (2004) for the region where the H 2 CO 
abundance is heavily enhanced because of evaporation of 
the grain mantles in low-mass protostars. This suggests a 
similar origin for the D2CO molecules in the hot core. 

Our previous single dish spectrum was too noisy to de- 
tect this wide emission (Paper II) . Recently, we have taken 
a good signal-to-noise ratio spectrum to compare with the 
interferometric observations (see Fig. 3). The single-dish 
spectrum shows two components, a narrow one at the ve- 
locity of the ambient cloud (already detected in Paper II) 
which is completely missed in our interferometric obser- 
vations and a wide one with a lincwidth similar to that of 
the PdB spectrum . Even in this wide component the in- 
terferometer recovers only the 50% of the flux, suggesting 
that part of them arises in a more extended area. 

A narrow and wide components with lincwidths similar 
to those of D 2 CO are also detected in the single dish spec- 
tra of the H 2 CO lines (Paper II). In Paper II, we have sep- 
arately mapped the narrow and wide components. These 
maps revealed that both components are associated with 
the outflow. While the wide component arises in the jet, 
the emission of the narrow component was more intense 
in the interface between the jet and the molecular cloud 
suggesting that it is tracing the molecular gas that is be- 
ing swept up by the jet. We propose that the missed flux 
in our interferometric D2CO image is probably arising in 
a more extended component associated with the bipolar 
outflow similarly to the case of H2CO. 

Our interferometric image shows that, in addition to 
the extended component, there is an intense and compact 
D2CO component associated with the hot core. We have 
estimated the D2CO abundance in this component assum- 
ing T rot =50 K and obtained N(D 2 CO)=3.5 10 14 cur 2 
and X(D2CO)~T.4 10~ 10 . This fractional abundance is a 
factor of ^2 higher than that derived in Paper II for the 
whole envelope. 



4.4. N 2 D + 

We have not detected the N2D + 3^2 line in our inter- 
ferometric image. The rms of the image is 50 mJy/beam 
which is about ~2% of the flux measured with the 30m 
telescope (Paper II). Thus, almost the 100% of the N 2 D+ 
emission has been missed in the interferometric observa- 
tions. This implies that the emission of N 2 D+ presents 
a quite uniform distribution across the envelope and sug- 
gests that it mainly arises in the cooler extended envelope. 
This is also in agreement with our interpretation in Paper 
II based on kinematical arguments. The linewidth of the 
N2D + 3^2 line, ~1 km s _1 , is a at least a factor of 2-3 
smaller than those of the lines arising in the warm com- 
ponent, and similar to those of the molecules which are 



Table 3. USB - 2.4 MHz 



Frcq. 


Molecule 


Transition 


Int." 


E„(K) 


231060.98 


ocs 


19 




18 


-2.6263 


100 


231187.69 


CH30CHO-E 


219,13 


— * 


218,14 


-4.8522 


179 


231199.96 


CH30CHO-A 


219,12 




218,13 


-4.8484 


179 


231200.14 


CH30CHO-E 


219,12 




218,13 


-4.8522 


179 


231220.99 


13 cs 


5 


— * 


4 


-4.2103 


51 


999 96 


N02 


4l,3 




4(1,4 


-4 8527 


12 


231239.07 


CH30CHO-A 


21g,i3 




218,14 


-4.8482 


179 


231232.06 




296,24 


- * 


295,25 


-5.0160 


253 


9^1 948 SI 




29 6 ,24 




294,26 


-5 3013 


253 


231252.98 


N02 


4l,3 


- * 


4o,4 


-5 0896 


12 


231 279 1 8 


CHqCHO-A 


12 6 ,7 


— * 


He, e 


-2.9783 


142 


231279.18 


CH3CHO-A 


12 6 ,6 


— * 


He, 5 


-2.9783 


142 


231286.40 


CH3OH 


10 2 ,2 




• 9 3 ,2 


-4.0689 


154 


231287.14 


NO2 


4l,3 


- * 


4o,4 


-5.0875 


12 


231310.42 


C2H5CN 


26i,25 


— * 


25i, 24 


-2.6517 


142 


231320.82 


CH3CHO-A 


125,8 




115,7 


-2.9001 


117 


231320.97 


CH3CHO-A 


125,7 


— * 


115,6 


-2.9001 


117 


231365.80 


C-C3D 


53,3 




43,2 


-3.1276 


28 


231410.27 


D 2 CO 


4 ,4 




4o,3 


-2.4242 


17 


231446.49 


CH3CHO-A 


124,9 




114,8 


-2.8386 


97 


231457.25 


CH3CHO-A 


124,8 




114,7 


-2.8386 


97 


231485.11 


c-C 3 D 


53,3 




4 3 .2 


-3.4082 


28 


231505.70 


HCOOH 


10i,g 




' 9l,8 


-3.1963 


53 


231511.81 


CH3OCHO-E 


274,24 




272,25 


-5.3849 


211 


231513.06 


CH3OCHO-E 


274,24 




273,25 


-5.1294 


211 


231558.55 


C2H5OH 


215,17 




214,18 


-3.6252 


215 


231559.62 


CH3OCHO-E 


27 5 ,23 




272,25 


-5.1292 


211 


231560.90 


C 2 H 5 OH 


20 5 ,i6 




20 4 ,i7 


-3.6250 


197 



a Base 10 logarithm of the integrated intensity in units of nm -MHz as 
appears tabulated in the JPL molecular line catalogue (Pickett et al. 



Table 4. LSB - 2.4 MHz 



Frcq. 


Molecule 


Transition 


Int" 


E„(K) 


228122.07 


CH3OCHO-A 


54,2 




4 2 ,3 


-7.1993 


8 


228134.42 


C-C3HD 


157,8 




157,9 


-3.9965 


288 


228147.97 


CH3OCHO-E 


18 6 ,13 




182,16 


-5.9576 


114 


228165.01 


CH3OCHO-A 


186,13 




184,14 


-5.9537 


114 


228165.86 


CH3OCHO-E 


44i4, 3 i 




43n,32 


-5.9853 


673 


228270.42 


CH3OCHO-A 


24 9 ,i6 




24 8 ,i7 


-4.8334 


220 


228272.26 


s 18 o 


9s 




8 8 


-3.9573 


82 


228288.69 


C2H5OH 


lln.o 




10 2 , 8 


-4.0647 


107 


228205.84 


CH3OCHO-E 


24 9 ,i6 




24 8 , 17 


-4.9577 


220 


228232 


Unidentified 












228245 


Unidentified 












228278.96 


HCCCHO 


25l,25 




24i, 24 


-2.8001 


135 


228310.56 


C2H3CHO 


132,12 




12i,n 


-4.6934 


36 


228358.20 


CH3OCHO-E 


24 9 ,i5 




24 8 , 17 


-5.4380 


220 


228415.69 


C-C3HD 


118, 4 




He, 5 


-4.5197 


169 


228434 


Unidentified 












228467 


Unidentified 












228483.13 


C2H5CN 


252,23 




242,22 


-5.4674 


158 


228494.88 


CH3OCHO-E 


144,10 




14l,l 3 


-6.2908 3 


62 


228522.71 


C2H5OH 


73,4 




72,6 


-4.3308 


84 


228544.17 


HCOOH 


102.8 




• 92,7 


-3.2324 


62 



a The same as in Table 3. 



thought to be good tracers of the cold envelope such as 
H 13 CO+ and N 2 H+. 



5. Chemical complexity of the hot core in 
NGC 7129-FIRS 2 

In Fig. 5 we show the observed wide band spectra in the 
upper side band (USB) and lower side band (LSB) of the 
1mm receiver. We have detected a real forest of lines in 
the 1mm band. The identification of these lines is not an 
easy task. Since we have moderate spectral resolution (2.5 
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USB (231.321 GHz) 



LSB (228.278 GHz) 



0.4 



>> 


m 
O 


I Q 




go" 


§ 0.2 






Li_ 


■f 


T I 



0.2 




231 .2 



228.5 



228.4 



228.3 



228.2 




231.4 



231 .2 



228.5 



228.4 



228.3 



228.2 



0.2 





231.4 



231 .2 



228.5 228.4 228.3 228.2 



0.2 




231.4 231.2 

Freq (GHz) 



228.5 



228.4 228.3 
Freq (GHz) 



228.2 



Fig. 5. Spectra of the upper side band (USB) and lower side band (LSB) of the 1mm receiver obtained in our PdBI 
observations with our best guess for line identifications. For clarity, we have divided the molecular species in three 
groups. The first group is formed by the more reliable identifications and S 18 0. Beginning from the top, in the first line 
of panels we compare the observed spectra with the synthesized ones taking into account only this first group. In the 
second line, we add CH 3 OCHO-A /E and show the obtained synthesized spectra. The inclusion of these compounds 
improve the agreement with the observational data. However, the large CH3OCHO-A /E abundance we derive from 
these observations, rise some questions about this identification. In the third line of panels, we show the synthesized 
spectra after adding some other exotic compounds. In the last line, it is shown the residual spectra after subtracting 
our fit. Note that the agreement between the synthesized and observed one is very good in the USB but some lines 
remain unidentified in the LSB. 
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Table 5. Fractional Abundances in the hot core of NGC 7129-FIRS 2 and comparison with other hot cores 





NGC 7129-FIRS 2 




NGC 1333 IRAS4 A 


IRAS 16293 


OMC 1 


G327.3-0.6 






~ 500 L 






~ 14 L 


- 23 L 


~1.0 10 4 L 


~ 10 5 L 


Molecule 


T * (K) 

±rot V / 


N*(cm~ 2 ) 


X a 




X" 


X c 


X d 


X e 


H 2 




7.9 xlO 24 


I 




1 


1 


1 


1 


CH 3 CN 


54 


5.5xl0 16 


7. Ox 10" 


-9 


1.6xl0~ 9 


1.0xl0~ 8 


4xl0~ 9 


7xl0" 7 


CH 3 OH 


50" 


2.4xl0 17 


3. Ox 10" 


-8 


<7 10" 9 


3.0xl0~ 7e 


lxlO" 7 


2 xlO" 5 


H 2 CO 


50" 


2.5 xlO 15 


3.1x10" 


10 


2.0xl0" 8 


6.0xl0" 8 


7.0xl0" 9 




D 2 CO 


50" 


3.5xl0 14 


4.4x10" 


11 




3.0xl0" 9s 






13 cs 


50 a 


3.4xl0 14 


4.3x10" 


11 






7xl0" n 


5xl0" 119 


ocs 


50 a 


2.4xl0 16 


3.0x10" 


-9 




l.OxlO -6 


7xl0" 8 


2xl0" 89 


HCOOH 


48 


3.5xl0 15 


4.4x10" 


10 


4.6 xlO" 9 


6.2xl0" 8 


8xl0" 10 


3X10" 1 '" 1 


C2H5OH 


75 


2.0xl0 16 


2.5x10" 


-9 






7xl0" 10 


3xl0" 9h 


C2H5CN 


50 a 


3.2xl0 15 


4.0x10" 


10 


<1.2xl0~ 9 


1.2xl0" 8 


3.0xl0" 9 


4xl0" 7 


C-C3D 


50 a 


2.5xl0 15 


3.1x10" 


10 










C-C3HD 


50 


1.8xl0 17 


2.3x10" 


-8 
















Uncertain identifications 








s 18 o 


50 a 


2.6xl0 16 


3.3x10" 


-9 




2.6xl0" 9 


3xl0" 10 


>6xl0 -13 


CH30CHO-E 


36 


4.0 xlO 18 


5.0x10" 


-7 


3.6xl0~ 8 


2.3xl0" 7 


lxlO" 8 


2xl0" 6 


CH30CHO-A 


36 a 


3.3xl0 ls 


4.1x10" 


-7 


3.4xl0" 8 


1.7xl0" 7 


lxlO" 8 


2xl0" 6 


CH3CHO-A 


50 a 


2.3xl0 15 


2.9x10" 


10 




1.9xl0" 8 




5xl0" 1()9 


HCCCHO 


50 a 


2.0 xlO 15 


2.5x10" 


10 










C2H3CHO 


50 a 


5.0xl0 16 


6.3x10" 


-9 










CH3OCH3 


50 a 


1.9xl0 15 


2.4x10" 


10 


<2.8xl0" s 


2.4xl0" 7 


8xl0" 9 


3xl0" 89 


34 so 2 


50 a 


3.8xl0 13 


4.8x10" 


12 




2.4xl0" 8 


2xl0" 10 


4xl0" 99 


no 2 


50 a 


l.OxlO 17 


1.2x10" 


-8 











* Values derived for the NGC 7129-FIRS 2 hot core assuming a size of ~0.72" x0.52" . To derive the H 2 column density we 
have considered that the hot core mass (2 M ) is uniformly distributed in an area ~a6 with a=0.004 pc (0.72") and fe=0.003 
pc (0.52"). 

a Assumed temperature. 

6 Values from Bottinelli et al. (2004) for a 0.0016 pc (0.5") source. 

c Values from Cazaux et al. (2003) and Wakelam et al. (2004) for a 0.0016 pc (2") source. 
d Values from Sutton et al. (1995) for a 0.034 pc (14") source. 
e Values from Gibb et al. (2000) for a 0.03 pc (2") source. 
f From Schoier et al. (2002). 
9 From Ceccarelli et al. (2001). 

h From Gibb et al. (2000) for a 0.32 pc (20") source. 



MHz), several candidates can be found for each line. On 
the other hand, in most cases we have two or more lines 
overlapped. We have tried to identify most of them with 
the following procedure: In a first step we select our best 
candidates for the lines and estimate the column density 
of the carriers assuming optically thin emission and lo- 
cal thermodynamic equilibrium (LTE). In a few cases, we 
are able to estimate the rotation temperature because we 
have at least two "isolated" lines (see Table 3). In the 
other cases, we have assumed T ro4 =50 K. In a second step, 
we synthesized the expected spectrum assuming the esti- 
mated column densities. The velocity profile is assumed to 
be Gaussian with a half power full width of 7 km s _1 and 
centered at the ambient cloud velocity, -10 km s _1 . Finally, 
we compare it with the observed spectrum and readjust 
the molecular column densities to improve the fit. We list 
in Tables 3 and 4 the lines detected in the USB and LSB 
spectra with the most likely line identification. 



Obviously, this identification procedure entails some 
ambiguity. We have classified the species according with 
the reliability of their identification in two groups that 
are differentiated in Table 5. The first group contains 
the most reliable identifications. In these cases, the ob- 
served frequencies are in total agreement with those cal- 
culated for these species. Moreover, the synthesized spec- 
trum matches quite well with the observed one (see top 
panel of Fig. 5). Finally, the derived fractional abundances 
are within the range of abundances measured for these 
species in hot cores. Within this group, we have OCS, 
13 CS, C 2 H 5 CN, c-C 3 D, D 2 CO, HCOOH, and C 2 H 5 OH in 
the USB and C 2 H 5 OH, C 3 HD in the LSB. The lines of 
the first group species are indicated by bold characters in 
Tables 3 and 4. The derived column densities and frac- 
tional abundances are shown in Table 5. 

We have put the complex molecules CH3OCHO-A 
and CH3OCHO-E (methyl formate) in the second group. 
Several lines of these species have been detected in the 
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USB and LSB (see second panels in Fig. 5). However, 
the derived fractional abundances in NGC 7129-FIRS 2 
are unexpectedly high compared to those measured in hot 
corinos and massive hot cores (see Table 3). One possi- 
bility is that the rotation temperature is higher than that 
assumed. To estimate the uncertainty due to the poorly 
known rotation temperature, we have repeated the column 
density calculations with T roi =300 K. The derived abun- 
dance is an order of magnitude lower and then it agrees 
with that derived in IRAS 16293. This high rotation tem- 
perature could be explained if this complex molecule arises 
in an inner and hotter region than CTI3CN. If methyl for- 
mate cannot be produced in the gas phase (as suggested 
by the laboratory results of Horn et al. 2004), our findings 
may imply that this species is formed on grain surfaces 
and that its binding energy is larger than that of CH 3 CN, 
thus requiring larger dust temperatures to desorb. 

The rarer isotope S 18 has also been included in the 
group of uncertain identifications. We have only one line 
of this molecule, S 18 9 8 ^8 8 at 228.272 GHz, and the de- 
rived S 18 abundance is unexpectedly high. In this case 
we could have a misidentification. There is a CH3OCHO-A 
line at 228.270 GHz, that could be an alternative identi- 
fication of the tentative S 18 line if the methyl formate 
abundance is unexpectedly high. 

The most uncertain identifications are C2H3CHO, 
CH3CHO-A, HCCCHO, and N0 2 . The observed frequen- 
cies are in agreement with those of C2H3CHO, CH3CHO- 
A, and HCCCHO lines if the emission arises in gas at 
a velocity v/ sr «— 5 km s _1 . This velocity is coincident 
with that of the bullet Rl (see Paper I) but is shifted 
by ^5 km s" 1 from that of the molecular cloud. Similarly, 
the fit to the USB spectrum improves if we include a high 
excitation methanol line (see Table 3) emitting at the ve- 
locities —5 and —20 km s _1 . These velocities are within 
the range of velocities associated with this Class proto- 
star but are shifted from that of the ambient cloud. NO2 
is not a common molecule in hot cores, and the three lines 
observed in the spectrum are overlapped with others of 
more common hot core species. We have kept it in the ta- 
ble because its inclusion improves the fit to the observed 
spectrum. 

In Fig. 5 we show the synthesized USB and LSB spec- 
tra assuming the molecular column densities listed in 
Table 3. The fit is quite good for the USB. However, 
it is poor for the LSB where we still have 4 unidenti- 
fied lines. This is quite frustrating taking into account 
that they are the most intense lines in this subband. 
The frequencies of the unidentified lines are 228467(±5), 
228434(±5),228245(±5) and 228232(±5) MHz, and we do 
not have any good candidates for them. There is a line of 
vibrationally excited DC3N at the frequency of 228467.44 
GHz. Emission of vibrationally excited HC3N has been de- 
tected in massive hot cores (see e.g. Wyrowski et al. 1999). 
However, if this identification were correct, we would ex- 
pect to see a vibrationally excited HC3N line at least ten 
times more intense at a frequency of 228303 MHz which 
is not detected. For the line at 228434 MHz our best can- 



didate is glycoaldehyde (CH 2 (OH)CHO). Glycoaldehyde 
has already been observed in the interstellar medium 
l(Hollis et al. 2004^ . However, the upper state energy of 
the line at 228434 MHz is too high (>1400 K) to be de- 
tectable in this source. The line at 228245 MHz could be 
propadienone (CH2CCO), but in this case we should have 
detected other lines of this compound in the same spec- 
trum. Finally, our best candidate for the line at 228232 
MHz is cyclopropenone (C3H 2 0). This molecule has never 
been detected in space. Since we have only one observable 
line in our spectra, we cannot confirm this detection. 

In addition to the lines detected in the wide-band spec- 
trum shown in Fig. 5, in Table 5 we also include the 
column densities derived for CH3OCH3 and 34 S02, these 
species detected in the CH3OH spectrum shown in Fig. 3. 
Since we have only one line for these species, the derived 
column density is also uncertain. 

6. Discussion 

In Table 5 we list the derived abundances of the molecules 
detected in the hot core towards FIRS 2. A first look 
at the Table reveals that FIRS 2 is rich in deuterated 
species (D2CO,c-C3D and C-C3HD), sulphuretted com- 
pounds ( 13 CS, OCS) and complex O- and N- molecules 
(HCOOH, C2H5OH, C2H5CN). Taking into account the 
large uncertainty in the estimated fractional abundances, 
we have centered our discussion in the most reliable iden- 
tifications, i.e., the first group in Table 3. 



6.1. Deuterated species 

The deuterium fractionation of molecular species is a well 
known process in dense interstellar clouds (Watson, W.D. 
1974). In the gas phase, this is mainly driven by the en- 
hancement of the H2D + /H 3 f ratio due to the exothermic 
proton-deuteron exchange reaction H^ + HD — > H2D 4 " 
+ H2 (+ 230 K), which does not proceed backwards if 
the gas temperature is less than about 30 K. This pro- 
cess is responsible for the deuteration of species such as 
HCO + and N2H + , via the reaction of H2D 4 " with the par- 
ent species CO and N?, respect i vely, as observed in d ense 
cloud cores l|Butner et al. 19951 IWilliams et al. 199S|) . 

Another phenomenon which boosts the deuterium 
fractionation is the freeze out of neutral species, in 
particular the abundant O and CO, which are effi- 



cient destruction partners of H3 and H2D H 



the 



gas phase. Thei r freeze out, known to happen in 
dense cold clouds Williams et al. 1995lCaselli et al.1999: 



iKramer et al. 19991: iBergin et al.2002HTafalla et al. 2002|) . 

further enhances the H2D + /H 3 f abundance ratio (be- 
cause of the enhanced formation rate of H2D + , and 
the reduced destruction rates for H2D+ and H^") 
and the conse quent transfer of t he deuteron to gas- 
phase species (ICaselli et al. 20021: iBacmann et al. 20031: 
ICrapsi et al. 2005|) . Once CO is highly depleted, mul- 
tiple deuterated forms of Hj~ can easily be produced 
(Roberts et al. 2003UWa,lmslev et al. 20oi. Indeed, large 
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abundances of H2D+ and D2H+ have been mea- 
sured toward two pre-stellar cores l|Caselli et al. 20o3 
IVastel et al. 2004^) . This has the consequence of en- 
hancing the gaseous D/H ratio to values around 
0.1, four orders of magnitude larger than the cos- 
mic D abundance, which allows an efficient deutera- 
tion on grain surfaces llTielens 19831 ICharnlev et al. 1997t 
IStantcheva fc Herbst 2003j) . 

We have detected the deuterated species D2CO, c- 
C 3 D and c-C 3 HD towards FIRS 2. Moreover, the D 2 CO 
interferometric image shows that the D2CO emission is 
strongly peaked towards the hot core position. Recent the- 
oretical and observational research suggests that the abun- 
dance of the non-deuterated species H 2 CO is enhanced by 
more than 2 orders of magnitude in the inner region of the 
envelope, where th e dust temperature is larger than 100 K 
l|Maret et al. 2004|) . The size of this region (-100 AU) is 
similar to that derived for the D2CO emission from the 
linewidth of the D2CO line. Assuming that evaporation is 
also the origin of D 2 CO, the [D 2 CO]/[H 2 CO] would tell us 
about the gaseous D/H ratio during the pre-stellar phase 
of the parent clump. Unfortunately, we have not interfer- 
ometric observations of H2CO and cannot derive an accu- 
rate [D 2 CO]/[H 2 CO] ratio in the hot core. Trying to get 
an estimate of the deuteration degree of formaldehyde, 
we have compared the single dish spectra of the D 2 CO 
4o4— >3o3 and H2CO 3i2— >2n lines (see Fig. 6). Both lines 
have similar excitation conditions and frequency. The line 
intensity ratio in the narrow component is lower by a 
factor of 2 than in the wide component suggesting a 
higher [T>2CO]/[H2CO] ratio in the hot core. Assuming 
the same physical conditions for D2CO and H2CO we 
derive [D 2 CO]/[H 2 CO]~0.14±0.03 in the hot core. This 
value is similar t o that found in pre-stellar cores and low- 
mass protostars dCaselli et al. 2002UBacmann et al. 2003; 
ICeccarelli et al. 200ll lLoinard et al. 2003h and it implies 
large D/H ratios during the cold pre-stellar phase 
(D/H~0.5, see Caselli et al. 2002b; Stantcheva & Herbst 
2003). 

Contrary to D2CO, we have not detected N2D + in our 
interferometric observations. It is interesting to note that 
whereas the deuteration of N2H" 1 " is a gas-phase process 
(which proceeds from the reaction between N2 and H2D + ), 
the formation of D2CO seems to require surface chemistry 
(|Ceccareili et al. 2002|) . Therefore, our result is not sur- 
prising. In fact, as stated above, the deuteration in the 
gas phase is efficient only if the gas temperature is lower 
than, say, 30 K. Thus, in a hot core, we do not expect to 
see N2D + . Moreover, in hot cores, where icy grain mantles 
evaporates, CO and H2O molecules are abundant in the 
gas phase and they efficiently destroy N 2 H + (and N 2 D+), 
so that one does not expect to see large abundances of 
N2H+ cither. The non-detection of N2D+ in the hot core 
associated with FIRS 2 confirms the interpretation given 
in Paper II that N2D 4 " mainly arises in the cool envelope, 
and constitutes an indirect proof of this chemical scheme. 
On the other hand, the large D2CO/H2CO abundance ra- 
tio observed in the same position can be easily explained 
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Fig. 6. Single-dish spectra of the H2CO 3i2— »2n and 
D2CO 4o4^3o3 lines towards FIRS 2. In the bottom panel 
we show the I(D 2 CO 4 4™>3o3)/I(H 2 CO 3 12 ->2 U ) line ra- 
tio. Note that this ratio presents a minimum around the 
velocity of the ambient cloud (~— 10 km s _1 ). 



if the deuterium fractionation of formaldehyde happened 
on grain surfaces during the cold pre-stellar phase, and 
if icy mantles have recently evaporated, so that gas phase 
chemistry did not have the time to significantly reduce the 
abundance of the deuterated "parent" species. 

6.2. Sulphuretted species 

The abundances of sulphuretted molecules are very depen- 
dent on the stellar age. We have compared the observed 
[CS]/[OCS] ratio in FIRS 2 (-1.1 assuming CS/ 13 CS=89) 
with the Viti et al. (2004) model for a 5 M star. The mea- 
sured ratio is consistent with the FIRS 2 being an evolved 
hot core with an age of >1.0 10 5 yr. The detection of large 
species like C2H5OH also suggests that we are dealing with 
an evolved hot core since these species are only formed at 
late times. Nevertheless, we would like to point out that 
the determination of the hot core age using the sulphuret- 
ted molecules is quite uncertain since other factors like 
the ice mantle composition and the hot core kinetic tem- 
perature can also strongly affect the abundances of these 
species (Wakelam et al. 2004). 
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6.3. Complex 0- and N-bearing molecules 

Theoretical models predict that the relative abundance 
of the N-bearing and O-bearing comp lex molecules is de- 
pendent on the hot core te mperature l|Caselli et al. 1 993: 
iRodeers fc Charnlev 200lh . The luminosity of FIRS 2 is 
a factor >10 larger than those of the two low mass hot 
corinos detected so far and a factor ^200 lower than that 
of the massive protostar creating the Orion hot core. The 
intermediate kinetic temperature of the FIRS 2 hot core 
is expected to produce a differentiated chemistry of com- 
plex molecules. In order to have a deeper insight into the 
dependence of the complex molecules abundances on the 
protostellar luminosity, in Table 5 we compare the ob- 
served abundances towards FIRS 2 with those measured 
in prototypical hot corinos and massive hot cores. The lu- 
minosity of the objects listed in Table 5 varies by 4 orders 
of magnitude, thus one expects to detect some chemical 
differences between these hot cores. 

However, the comparison of abundances between dif- 
ferent hot cores (corinos) is difficult. First of all, since mas- 
sive stars are usually located farther from the Sun than 
their low- and intermediate-mass counterparts, the abun- 
dances are derived in different spatial scales. The frac- 
tional abundances shown in Table 5 are averaged values 
in regions with sizes ranging from ~0.002 pc in NGC1333- 
IRAS 4A to 0.32 pc in OMC1 and G327. These sizes are 
given by the angular resolution of the current instrumen- 
tation at the hot core distance. Thus, the derived abun- 
dances in massive hot cores would be severely underes- 
timated in the case of sizes similar to those of hot cori- 
nos. On the other hand, we can have several hot cores in 
the studied region, specially in massive star forming re- 
gions where clustering is more important and are located 
at greater distances. In order to minimize these problems, 
we normalize the studied abundances to those of the par- 
ent species. 

As a first step, we have normalized all the molecular 
abundances to that of H2CO and CH3OH. These species 
are thought to be the "parent" molecules of complex O- 
bearing species, although some do ubts are cast in the 
particular case of methyl formate l|Horn et al. 2004!) . In 
Fig. 7a, we show the relative abundances of HCOOH, 
CH3OH, CH 3 CN and C2H5CN respect to H 2 CO as a func- 
tion of the luminosity. Since we have not interferomet- 
ric observations of H2CO in FIRS 2, we have estimated 
the H 2 CO abundance assuming [D 2 CO]/[H 2 CO]=0.014. 
A clear trend is observed in the relative abundance of 
CH 3 OH, CH 3 CN and C 2 H 5 CN respect to H 2 CO as a func- 
tion of the stellar luminosity (Fig. 7a). The abundance 
of all these molecules increases by a factor of ^10 from 
NGC 1333-IRAS 4A (L=14 L Q ) to OMC1 (L=10 4 L ). 
This trend does not present significant differences between 
O- and N-bearing molecules. Contrary to these molecules, 
the [HCOOH] /[H 2 CO] ratio does not present a systematic 
behavior. 

In Fig. 7b we plot the abundances of the same 
molecules normalized to that of CH3OH. Different be- 



haviors are found for the different molecules. The 
[CH 3 CN]/[CH 3 OH] ratio does not present any systematic 
trend with the protostellar luminosity, although the dis- 
persion in the values is quite high. The same remains true 
for C 2 H 5 CN and probably for C 2 H 5 OH (we only have 3 
points for C 2 H5 0H). This suggests that the ratio between 
O- and N-bearing complex molecules does not seem to 
strongly depend on the protostellar luminosity, contrary to 
that expected from theoretical models. However, the rela- 
tive abundance of HCOOH seems to decrease by 2 orders 
of magnitude from NGC 1333-IRAS4 A to G327.3. Thus, 
the relative abundances of HCOOH to those of CH3OH, 
CH 3 CN, C 2 H 5 CN and C 2 H 5 OH seem to decrease with the 
protostellar luminosity. 

We propose that the increase of CH3 0H/H 2 CO 
(as well as CH3OH/HCOOH), CH 3 CN/H 2 CO, and 
C 2 H5CN/H 2 CO with stellar luminosity is likely due to 
differences in the grain mantle composition between low 
and massive star forming regions. Indeed, differences in 
the initial conditions (in particular gas density and dust 
temperature) are known to strongly affect surface chem- 
istry. For example, Caselli et al. (1993) found that complex 
N-bearing species are more easily formed if the dust grain 
temperature (Tdust) is about 40 K, during the collapse 
phase. Although such temperature can be found in a sig- 
nificant fraction of a collapsing massive clump surrou nding 
a massive young stellar object ijFontani et al. 20*02^) . this 
is not the case for low mass cores. Therefore, the observed 
trend (at least for N-bearing species) is consistent with 
theoretical predictions. But Caselli et al. (1993) also found 
that at Trf„ st ~40 K, methanol is not efficiently formed 
on the surface because it requires the volatile H to stick 
on the grain surface (see also Charnley et al. 1992; van 
der Tak et al. 2000). Thus, the observed CH 3 OH/H 2 CO 
increase with stellar luminosity is not well understood, 
unless another surface formation mechanism for methanol 
(besides the CO hydrogenation) becomes important in the 
warmer high mass clumps, e.g. OH+CH3 (see also discus- 
sion in Pontoppidan et al. 2003). However, this conclusion 
is only based on 4 sources and has to be confirmed with a 
larger and more complete sample. Observational bias such 
as the different spatial scale for the different objects can 
contribute to give this trend if the spatial extension of the 
emission of the observed molecules is different. 

In Fig. 7c, we compare only the N-bearing molecules. 
The [C 2 H5CN]/ [CH3CN] ratio remains quite constant 
with a dispersion of about a factor <10 between all 
the sources. This uniform [C 2 H 5 CN]/ [CH 3 CN] ratio 
suggesting that the chemistry of both compounds is 
linked. Indeed, these two molecules have a common parent 
species: C 2 N (see e.g. Fig. 7 in Caselli et al. 1993) 

7. Conclusions: The chemistry of the IM hot core 
NGC 7129-FIRS 2 

In this paper we present the first detection of an 
intermediate-mass hot core. Our intcrferometric obser- 
vations towards FIRS 2 allow us to estimate a size of 



A. Fuente et al.: Detection of a hot core in the intermediate-mass Class protostar NGC 7129-FIRS 2 



13 



650x900 AU and a mass of 2 M Q for the hot core as- 
sociated with this IM object. The dimensions and mass 
of this IM hot core are intermediate between those mea- 
sured in hot corinos (r^l50 AU,M<1 M Q ) and massive 
stars (r^3000 AU,M>10 M Q ) and consequently one ex- 
pects to find a differentiated chemistry in it. 

Our interferometric observations of CH3CN, CH3OH, 
D2CO and N2D" 1 " provides information on the structure 
and chemistry of this object. While the CH3CN emis- 
sion arises mainly in the hot core, the CH3OH emission 
has a component associated with the bipolar outflow. The 
CH3CN abundance is ~ 7.0 10~ 9 in the hot core , which 
is about 3 orders of magnitude larger than in the cool en- 
velope. The methanol abundance is enhanced (>3 10 _8 -a 
few 10~ 7 ) in the hot core and the outflow components. 
However the rotation temperature is higher in the hot 
core than in the outflow revealing very different physical 
conditions, and probably a different CH3OH desorption 
mechanism for the two components. 

The doubly deuterated formaldehyde also presents 
enhanced abundances in the hot core. We estimate a 
[D 2 CO]/[H 2 CO]~0.14, which is 4 orders of magnitude 
larger than the cosmic D abundance and similar to those 
found in pre-stellar clumps and low-mass protostars. The 
enhanced [D 2 CO]/[H 2 CO] abundance in the hot core sug- 
gests that grain surface chemistry is responsible of the 
deuteration process of H 2 CO. In contrast to the high 
deuteration degree of formaldehyde, we have not detected 
N 2 D + in the hot core. This is consistent with the chemical 
scheme in which while the deuteration of H 2 CO requires 
surface chemistry, the deuteration of N 2 H + is a gas-phase 
process. 

A large number of molecular lines have been de- 
tected in our interferometric spectra towards FIRS 2. 
Most of these lines are identified as belonging to 
deuterated (D 2 CO,c-C3D and C-C3HD), sulphuret- 
ted ( 13 CS,OCS) and complex 0-/N-bearing species 
(HCOOH,C 2 H 5 OH,C 2 H 5 CN). Deuterated species whose 
deuteration requires surface chemistry, such as D 2 CO, 
present enhanced abundances in the warm regions asso- 
ciated with low-mass protostars. Loinard et al. (2003) 
searched for the doubly deuterated form of formalde- 
hyde (D 2 CO) in a large sample of young stellar ob- 
jects. D 2 CO was detected in all low-mass protostars, with 
[D 2 CO]/[H 2 CO] ratios of 2-40%. On the other hand, no 
detection was obtained towards more massive protostars 
(where [D 2 CO] [H 2 CO]<0.5%). This is consistent with the 
results reported by Turner (1990) who detected D 2 CO in 
Orion and measured [D 2 CO]/[H 2 CO]~0.003. If the hot 
cores associated with massive stars are older or signifi- 
cantly denser than those surrounding low-mass objects, 
gas phase chemistry could have had the time to re-set the 
deuterium fractionation to values close to the cosmic D /H 
ratio. An alternative explanation is that the gas temper- 
ature of the mate rial accreting high mass protostars is 
larger than 30 K ijFontani et al. 2002]) . so that the deu- 
terium fractionation efficiency in the gas phase is strongly 
reduced already before the hot core phase 



The sulphuretted and complex compounds are charac- 
teristic of hot cores in both the low-mass and the high- 
mass regimes. FIRS 2 is the first IM hot core detected 
thus far and presents a unique opportunity to study the 
link between the chemistry of hot corinos and massive 
hot cores. We have compared the abundances of complex 
molecules in FIRS 2 with those in hot corinos and the mas- 
sive hot cores OMC1 and G327.3-0.6. Contrary to model 
predictions, we do no detect any dependence of the O- 
/N-complex molecules ratio on the protostellar luminos- 
ity. However, we detect differences between the behavior 
of the O-bearing species with the stellar luminosity. While 
H 2 CO and HCOOH are more abundant in low luminosity 
sources, CH3OH seems to be more abundant in massive 
objects. We propose that this could be due to a different 
mantle composition in the two kinds of region, caused by 
different physical conditions (mainly gas density and dust 
temperature) during the pre-stellar and accretion phase. 
However, this could also be due to other factors such as 
the different spatial scale of the observations or a possible 
contribution of the shocked gas associated with the bipo- 
lar outflow to the emission of these molecules. Differences 
in the hot core ages and/or cloud initial conditions could 
also produce these differences. The detection and detailed 
study of more intermediate-mass hot cores are necessary 
to establish firm conclusions. 
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